APOBEC shapes tumor evolution and age at onset of lung cancer in smokers.

APOBEC enzymes are part of the innate immunity and are responsible for restricting viruses and retroelements by deaminating cytosine residues1,2. Most solid tumors harbor different levels of somatic mutations attributed to the off-target activities of APOBEC3A (A3A) and/or APOBEC3B (A3B)3-6. However, how APOBEC3A/B enzymes shape the tumor evolution in the presence of exogenous mutagenic processes is largely unknown. Here, by combining deep whole-genome sequencing with multi-omics profiling of 309 lung cancers from smokers with detailed tobacco smoking information, we identify two subtypes defined by low (LAS) and high (HAS) APOBEC mutagenesis. LAS are enriched for A3B-like mutagenesis and KRAS mutations, whereas HAS for A3A-like mutagenesis and TP53 mutations. Unlike APOBEC3A, APOBEC3B expression is strongly associated with an upregulation of the base excision repair pathway. Hypermutation by unrepaired A3A and tobacco smoking mutagenesis combined with TP53-induced genomic instability can trigger senescence7, apoptosis8, and cell regeneration9, as indicated by high expression of pulmonary healing signaling pathway, stemness markers and distal cell-of-origin in HAS. The expected association of tobacco smoking variables (e.g., time to first cigarette) with genomic/epigenomic changes are not observed in HAS, a plausible consequence of frequent cell senescence or apoptosis. HAS have more neoantigens, slower clonal expansion, and older age at onset compared to LAS, particularly in heavy smokers, consistent with high proportions of newly generated, unmutated cells and frequent immuno-editing. These findings show how heterogeneity in mutational burden across co-occurring mutational processes and cell types contributes to tumor development, with important clinical implications.

Glycidamide-induced hypermutation in yeast single-stranded DNA reveals a ubiquitous clock-like mutational motif in humans.

Mutagens often prefer specific nucleotides or oligonucleotide motifs that can be revealed by studying the hypermutation spectra in single-stranded (ss) DNA. We utilized a yeast model to explore mutagenesis by glycidamide, a simple epoxide formed endogenously in humans from the environmental toxicant acrylamide. Glycidamide caused ssDNA hypermutation in yeast predominantly in cytosines and adenines. The most frequent mutations in adenines occurred in the nAt→nGt trinucleotide motif. Base substitutions A→G in this motif relied on Rev1 translesion polymerase activity. Inactivating Rev1 did not alter the nAt trinucleotide preference, suggesting it may be an intrinsic specificity of the chemical reaction between glycidamide and adenine in the ssDNA. We found this mutational motif enriched in published sequencing data from glycidamide-treated mouse cells and ubiquitous in human cancers. In cancers, this motif was positively correlated with the single base substitution (SBS) smoking-associated SBS4 signature, with the clock-like signatures SBS1, SBS5, and was strongly correlated with smoking history and with age of tumor donors. Clock-like feature of the motif was also revealed in cells of human skin and brain. Given its pervasiveness, we propose that this mutational motif reflects mutagenic lesions to adenines in ssDNA from a potentially broad range of endogenous and exogenous agents.

Changes in metabolic landscapes shape divergent but distinct mutational signatures and cytotoxic consequences of redox stress.

Mutational signatures discerned in cancer genomes, in aging tissues and in cells exposed to toxic agents, reflect complex processes underlying transformation of cells from normal to dysfunctional. Due to its ubiquitous and chronic nature, redox stress contributions to cellular makeover remain equivocal. The deciphering of a new mutational signature of an environmentally-relevant oxidizing agent, potassium bromate, in yeast single strand DNA uncovered a surprising heterogeneity in the mutational signatures of oxidizing agents. NMR-based analysis of molecular outcomes of redox stress revealed profound dissimilarities in metabolic landscapes following exposure to hydrogen peroxide versus potassium bromate. The predominance of G to T substitutions in the mutational spectra distinguished potassium bromate from hydrogen peroxide and paraquat and mirrored the observed metabolic changes. We attributed these changes to the generation of uncommon oxidizing species in a reaction with thiol-containing antioxidants; a nearly total depletion of intracellular glutathione and a paradoxical augmentation of potassium bromate mutagenicity and toxicity by antioxidants. Our study provides the framework for understanding multidimensional processes triggered by agents collectively known as oxidants. Detection of increased mutational loads associated with potassium bromate-related mutational motifs in human tumors may be clinically relevant as a biomarker of this distinct type of redox stress.

Genomic and evolutionary classification of lung cancer in never smokers.

Lung cancer in never smokers (LCINS) is a common cause of cancer mortality but its genomic landscape is poorly characterized. Here high-coverage whole-genome sequencing of 232 LCINS showed 3 subtypes defined by copy number aberrations. The dominant subtype (piano), which is rare in lung cancer in smokers, features somatic UBA1 mutations, germline AR variants and stem cell-like properties, including low mutational burden, high intratumor heterogeneity, long telomeres, frequent KRAS mutations and slow growth, as suggested by the occurrence of cancer drivers' progenitor cells many years before tumor diagnosis. The other subtypes are characterized by specific amplifications and EGFR mutations (mezzo-forte) and whole-genome doubling (forte). No strong tobacco smoking signatures were detected, even in cases with exposure to secondhand tobacco smoke. Genes within the receptor tyrosine kinase-Ras pathway had distinct impacts on survival; five genomic alterations independently doubled mortality. These findings create avenues for personalized treatment in LCINS.

UV-exposure, endogenous DNA damage, and DNA replication errors shape the spectra of genome changes in human skin.

Human skin is continuously exposed to environmental DNA damage leading to the accumulation of somatic mutations over the lifetime of an individual. Mutagenesis in human skin cells can be also caused by endogenous DNA damage and by DNA replication errors. The contributions of these processes to the somatic mutation load in the skin of healthy humans has so far not been accurately assessed because the low numbers of mutations from current sequencing methodologies preclude the distinction between sequencing errors and true somatic genome changes. In this work, we sequenced genomes of single cell-derived clonal lineages obtained from primary skin cells of a large cohort of healthy individuals across a wide range of ages. We report here the range of mutation load and a comprehensive view of the various somatic genome changes that accumulate in skin cells. We demonstrate that UV-induced base substitutions, insertions and deletions are prominent even in sun-shielded skin. In addition, we detect accumulation of mutations due to spontaneous deamination of methylated cytosines as well as insertions and deletions characteristic of DNA replication errors in these cells. The endogenously induced somatic mutations and indels also demonstrate a linear increase with age, while UV-induced mutation load is age-independent. Finally, we show that DNA replication stalling at common fragile sites are potent sources of gross chromosomal rearrangements in human cells. Thus, somatic mutations in skin of healthy individuals reflect the interplay of environmental and endogenous factors in facilitating genome instability and carcinogenesis.

Somatic mutation distributions in cancer genomes vary with three-dimensional chromatin structure.

Somatic mutations in driver genes may ultimately lead to the development of cancer. Understanding how somatic mutations accumulate in cancer genomes and the underlying factors that generate somatic mutations is therefore crucial for developing novel therapeutic strategies. To understand the interplay between spatial genome organization and specific mutational processes, we studied 3,000 tumor-normal-pair whole-genome datasets from 42 different human cancer types. Our analyses reveal that the change in somatic mutational load in cancer genomes is co-localized with topologically-associating-domain boundaries. Domain boundaries constitute a better proxy to track mutational load change than replication timing measurements. We show that different mutational processes lead to distinct somatic mutation distributions where certain processes generate mutations in active domains, and others generate mutations in inactive domains. Overall, the interplay between three-dimensional genome organization and active mutational processes has a substantial influence on the large-scale mutation-rate variations observed in human cancers.

Similarity between mutation spectra in hypermutated genomes of rubella virus and in SARS-CoV-2 genomes accumulated during the COVID-19 pandemic.

Genomes of tens of thousands of SARS-CoV2 isolates have been sequenced across the world and the total number of changes (predominantly single base substitutions) in these isolates exceeds ten thousand. We compared the mutational spectrum in the new SARS-CoV-2 mutation dataset with the previously published mutation spectrum in hypermutated genomes of rubella-another positive single stranded (ss) RNA virus. Each of the rubella virus isolates arose by accumulation of hundreds of mutations during propagation in a single subject, while SARS-CoV-2 mutation spectrum represents a collection events in multiple virus isolates from individuals across the world. We found a clear similarity between the spectra of single base substitutions in rubella and in SARS-CoV-2, with C to U as well as A to G and U to C being the most prominent in plus strand genomic RNA of each virus. Of those, U to C changes universally showed preference for loops versus stems in predicted RNA secondary structure. Similarly, to what was previously reported for rubella virus, C to U changes showed enrichment in the uCn motif, which suggested a subclass of APOBEC cytidine deaminase being a source of these substitutions. We also found enrichment of several other trinucleotide-centered mutation motifs only in SARS-CoV-2-likely indicative of a mutation process characteristic to this virus. Altogether, the results of this analysis suggest that the mutation mechanisms that lead to hypermutation of the rubella vaccine virus in a rare pathological condition may also operate in the background of the SARS-CoV-2 viruses currently propagating in the human population.

Mixture Analyses of Air-sampled Pollen Extracts Can Accurately Differentiate Pollen Taxa.

The daily pollen forecast provides crucial information for allergic patients to avoid exposure to specific pollen. Pollen counts are typically measured with air samplers and analyzed with microscopy by trained experts. In contrast, this study evaluated the effectiveness of identifying the component pollens using the metabolites extracted from an air-sampled pollen mixture. Ambient air-sampled pollen from Munich in 2016 and 2017 was visually identified from reference pollens and extracts were prepared. The extracts were lyophilized, rehydrated in optimal NMR buffers, and filtered to remove large proteins. NMR spectra were analyzed for pollen associated metabolites. Regression and decision-tree based algorithms using the concentration of metabolites, calculated from the NMR spectra outperformed algorithms using the NMR spectra themselves as input data for pollen identification. Categorical prediction algorithms trained for low, medium, high, and very high pollen count groups had accuracies of 74% for the tree, 82% for the grass, and 93% for the weed pollen count. Deep learning models using convolutional neural networks performed better than regression models using NMR spectral input, and were the overall best method in terms of relative error and classification accuracy (86% for tree, 89% for grass, and 93% for weed pollen count). This study demonstrates that NMR spectra of air-sampled pollen extracts can be used in an automated fashion to provide taxa and type-specific measures of the daily pollen count.

Similarity between mutation spectra in hypermutated genomes of rubella virus and in SARS-CoV-2 genomes accumulated during the COVID-19 pandemic.

Genomes of tens of thousands of SARS-CoV2 isolates have been sequenced across the world and the total number of changes (predominantly single base substitutions) in these isolates exceeds ten thousand. We compared the mutational spectrum in the new SARS-CoV-2 mutation dataset with the previously published mutation spectrum in hypermutated genomes of rubella - another positive single stranded (ss) RNA virus. Each of the rubella isolates arose by accumulation of hundreds of mutations during propagation in a single subject, while SARS-CoV-2 mutation spectrum represents a collection events in multiple virus isolates from individuals across the world. We found a clear similarity between the spectra of single base substitutions in rubella and in SARS-CoV-2, with C to U as well as A to G and U to C being the most prominent in plus strand genomic RNA of each virus. Of those, U to C changes universally showed preference for loops versus stems in predicted RNA secondary structure. Similarly, to what was previously reported for rubella, C to U changes showed enrichment in the uCn motif, which suggested a subclass of APOBEC cytidine deaminase being a source of these substitutions. We also found enrichment of several other trinucleotide-centered mutation motifs only in SARS-CoV-2 - likely indicative of a mutation process characteristic to this virus. Altogether, the results of this analysis suggest that the mutation mechanisms that lead to hypermutation of the rubella vaccine virus in a rare pathological condition may also operate in the background of the SARS-CoV-2 viruses currently propagating in the human population.

APOBEC3-dependent kataegis and TREX1-driven chromothripsis during telomere crisis.

Chromothripsis and kataegis are frequently observed in cancer and may arise from telomere crisis, a period of genome instability during tumorigenesis when depletion of the telomere reserve generates unstable dicentric chromosomes1-5. Here we examine the mechanism underlying chromothripsis and kataegis by using an in vitro telomere crisis model. We show that the cytoplasmic exonuclease TREX1, which promotes the resolution of dicentric chromosomes4, plays a prominent role in chromothriptic fragmentation. In the absence of TREX1, the genome alterations induced by telomere crisis primarily involve breakage-fusion-bridge cycles and simple genome rearrangements rather than chromothripsis. Furthermore, we show that the kataegis observed at chromothriptic breakpoints is the consequence of cytosine deamination by APOBEC3B. These data reveal that chromothripsis and kataegis arise from a combination of nucleolytic processing by TREX1 and cytosine editing by APOBEC3B.

Analysis of APOBEC-induced mutations in yeast strains with low levels of replicative DNA polymerases.

Yeast strains with low levels of the replicative DNA polymerases (alpha, delta, and epsilon) have high levels of chromosome deletions, duplications, and translocations. By examining the patterns of mutations induced in strains with low levels of DNA polymerase by the human protein APOBEC3B (a protein that deaminates cytosine in single-stranded DNA), we show dramatically elevated amounts of single-stranded DNA relative to a wild-type strain. During DNA replication, one strand (defined as the leading strand) is replicated processively by DNA polymerase epsilon and the other (the lagging strand) is replicated as short fragments initiated by DNA polymerase alpha and extended by DNA polymerase delta. In the low DNA polymerase alpha and delta strains, the APOBEC-induced mutations are concentrated on the lagging-strand template, whereas in the low DNA polymerase epsilon strain, mutations occur on the leading- and lagging-strand templates with similar frequencies. In addition, for most genes, the transcribed strand is mutagenized more frequently than the nontranscribed strand. Lastly, some of the APOBEC-induced clusters in strains with low levels of DNA polymerase alpha or delta are greater than 10 kb in length.

Correction: Repair of multiple simultaneous double-strand breaks causes bursts of genome-wide clustered hypermutation.

[This corrects the article DOI: 10.1371/journal.pbio.3000464.].

Mutation signatures specific to DNA alkylating agents in yeast and cancers.

Alkylation is one of the most ubiquitous forms of DNA lesions. However, the motif preferences and substrates for the activity of the major types of alkylating agents defined by their nucleophilic substitution reactions (SN1 and SN2) are still unclear. Utilizing yeast strains engineered for large-scale production of single-stranded DNA (ssDNA), we probed the substrate specificity, mutation spectra and signatures associated with DNA alkylating agents. We determined that SN1-type agents preferably mutagenize double-stranded DNA (dsDNA), and the mutation signature characteristic of the activity of SN1-type agents was conserved across yeast, mice and human cancers. Conversely, SN2-type agents preferably mutagenize ssDNA in yeast. Moreover, the spectra and signatures derived from yeast were detectable in lung cancers, head and neck cancers and tumors from patients exposed to SN2-type alkylating chemicals. The estimates of mutation loads associated with the SN2-type alkylation signature were higher in lung tumors from smokers than never-smokers, pointing toward the mutagenic activity of the SN2-type alkylating carcinogens in cigarettes. In summary, our analysis of mutations in yeast strains treated with alkylating agents, as well as in whole-exome and whole-genome-sequenced tumors identified signatures highly specific to alkylation mutagenesis and indicate the pervasive nature of alkylation-induced mutagenesis in cancers.

Comprehensive analysis of chromothripsis in 2,658 human cancers using whole-genome sequencing.

Chromothripsis is a mutational phenomenon characterized by massive, clustered genomic rearrangements that occurs in cancer and other diseases. Recent studies in selected cancer types have suggested that chromothripsis may be more common than initially inferred from low-resolution copy-number data. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we analyze patterns of chromothripsis across 2,658 tumors from 38 cancer types using whole-genome sequencing data. We find that chromothripsis events are pervasive across cancers, with a frequency of more than 50% in several cancer types. Whereas canonical chromothripsis profiles display oscillations between two copy-number states, a considerable fraction of events involve multiple chromosomes and additional structural alterations. In addition to non-homologous end joining, we detect signatures of replication-associated processes and templated insertions. Chromothripsis contributes to oncogene amplification and to inactivation of genes such as mismatch-repair-related genes. These findings show that chromothripsis is a major process that drives genome evolution in human cancer.

The repertoire of mutational signatures in human cancer.

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature1. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses3-15, enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated-but distinct-DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Infectious vaccine-derived rubella viruses emerge, persist, and evolve in cutaneous granulomas of children with primary immunodeficiencies.

Rubella viruses (RV) have been found in an association with granulomas in children with primary immune deficiencies (PID). Here, we report the recovery and characterization of infectious immunodeficiency-related vaccine-derived rubella viruses (iVDRV) from diagnostic skin biopsies of four patients. Sequence evolution within PID hosts was studied by comparison of the complete genomic sequences of the iVDRVs with the genome of the vaccine virus RA27/3. The degree of divergence of each iVDRV correlated with the duration of persistence indicating continuous intrahost evolution. The evolution rates for synonymous and nonsynonymous substitutions were estimated to be 5.7 x 10-3 subs/site/year and 8.9 x 10-4 subs/site/year, respectively. Mutational spectra and signatures indicated a major role for APOBEC cytidine deaminases and a secondary role for ADAR adenosine deaminases in generating diversity of iVDRVs. The distributions of mutations across the genes and 3D hotspots for amino acid substitutions in the E1 glycoprotein identified regions that may be under positive selective pressure. Quasispecies diversity was higher in granulomas than in recovered infectious iVDRVs. Growth properties of iVDRVs were assessed in WI-38 fibroblast cultures. None of the iVDRV isolates showed complete reversion to wild type phenotype but the replicative and persistence characteristics of iVDRVs were different from those of the RA27/3 vaccine strain, making predictions of iVDRV transmissibility and teratogenicity difficult. However, detection of iVDRV RNA in nasopharyngeal specimen and poor neutralization of some iVDRV strains by sera from vaccinated persons suggests possible public health risks associated with iVDRV carriers. Detection of IgM antibody to RV in sera of two out of three patients may be a marker of virus persistence, potentially useful for identifying patients with iVDRV before development of lesions. Studies of the evolutionary dynamics of iVDRV during persistence will contribute to development of infection control strategies and antiviral therapies.

Repair of multiple simultaneous double-strand breaks causes bursts of genome-wide clustered hypermutation.

A single cancer genome can harbor thousands of clustered mutations. Mutation signature analyses have revealed that the origin of clusters are lesions in long tracts of single-stranded (ss) DNA damaged by apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) cytidine deaminases, raising questions about molecular mechanisms that generate long ssDNA vulnerable to hypermutation. Here, we show that ssDNA intermediates formed during the repair of gamma-induced bursts of double-strand breaks (DSBs) in the presence of APOBEC3A in yeast lead to multiple APOBEC-induced clusters similar to cancer. We identified three independent pathways enabling cluster formation associated with repairing bursts of DSBs: 5' to 3' bidirectional resection, unidirectional resection, and break-induced replication (BIR). Analysis of millions of mutations in APOBEC-hypermutated cancer genomes revealed that cancer tolerance to formation of hypermutable ssDNA is similar to yeast and that the predominant pattern of clustered mutagenesis is the same as in resection-defective yeast, suggesting that cluster formation in cancers is driven by a BIR-like mechanism. The phenomenon of genome-wide burst of clustered mutagenesis revealed by our study can play an important role in generating somatic hypermutation in cancers as well as in noncancerous cells.